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Abstract 



We study lepton flavour-violating interactions which could result in the r-lepton 
£^ ' production in the v^N scattering or in \x — ► r conversion on nucleons at high ener- 

gies. Phenomenological bounds on the strength of fv^qq interactions are extracted 
from the combined result of the NOMAD and CHORUS experiments on searching 
for — > v T oscillations. Some of these bounds supersede limits from rare decays. 

We also propose a "missing energy" type experiment searching for fj, — ► r conver- 
sion on nucleons. The experiment can be performed at a present accelerator or at a 
^ ' future neutrino factory. 
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1 Introduction 



In the standard model of strong and electroweak interactions (SM) with mass- 
less neutrinos lepton flavour is conserved. However many extensions of the 
SM including the simplest extension with massive neutrinos predict flavour- 
violating interactions. For instance leptoquarks which occur naturally in all 
unified models such as SU(5) [1] or Pati-Salam 5 , f/(4) [2] , in models with 
quark-lepton substructure [3] and in supersymmetric models with R-parity 
violation [4]. In general leptoquarks have baryon or lepton number violating 
couplings and they must be very heavy in order to avoid rapid proton decay or 
large Majorana neutrino masses. However, in theories with conserved baryon 
and lepton number leptoquark masses and couplings satisfy much weaker 
bounds which come mainly from the analysis of flavour-violating processes. 
Many searches for specific reactions, namely fi — e conversion on nuclei [7], 
/U — > e7 [5] and fi — > 3e [6] decays, rare K [8], D- B- [9] and r— decays [10] 
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and also search for lepton flavour violation in e — p collisions [12] and neu- 
trino interactions [13] have been performed resulting in stringent limits on the 
coupling strength of lepton flavour violating four-fermion 11 qq interactions. 

There are also projects to improve some of these bounds significantly, see e.g. 
[14,15] including experiments at a neutrino factory. The neutrino factory is a 
future powerful tool for investigation of the nature of neutrino, in particular 
neutrino osccilations. It could also be useful for investigation of other physics 
beyond the standard model, see e.g. [18] and references therein. This facility 
would be an intense source of muons and might enable to improve the current 
upper limits on flavour-violating transitions, involving in particular /i's or z/^'s, 
by several orders of magnitude. 

In this note we have suggested to probe lepton flavour violation through 
r-lepton production in u^N scattering or p + N — > r + N conversion on 
nucleons at high energies. It should be stressed that experiments on search 
for p-r conversion will probe flavour-violating interactions, in particular like 
L = pr[G\uc + G2<lb + G^sb]. Due to trivial kinematical constraints (e.g. 
r — > p + D decay cannot occur, because of m T < m M + m D ), there are no 
stringent bounds on such interaction from the existing experimental bounds 
on rare r-lepton, such as r — > p + 7 [10], and charm/beauty meson decay 
rates [9]. For investigation of such flavour- violating interactions to perform 
experiments on searching for p-r conversion is important. 

This paper is organised as follows. In Sec. 2 we extract phenomenological 
bounds on the strength of fu^qq interactions from the current experimen- 
tal bounds on search for — v T oscillations from the NOMAD and CHORUS 
experiments. These bounds can compete from bounds derived from r-decays. 
In Sec. 3 we consider a "missing energy" type experiment to search for lepton 
flavour violation in p — r conversion on nucleons at high energies. In Sec. 4 
we present our conclusions. 



2 Phenomenological bounds 

The latest combined limit of the CHORUS and NOMAD experiments on the 
probability of — >• v T oscillations is [19]: 



PW ~* V r) < 0.6 • 10- 4 (1) 
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Using the bound of Eq.(l) one can find that NOMAD data lead to the bound 

on flavour- violating cross section of neutrino with isoscalar nucleon target. 

Bound of Eq.(2) allows to obtain restriction on flavour- violating interactions 
fv^qq which lead to nonzero cross section for the flavour changing reaction 
v^N — > t + .... Let us consider the case of scalar leptoquarks. SU(3)®SU(2)<g> 
U(l) invariant Lagrangian of scalar and vector leptoquarks satisfying baryon 
and lepton number conservation has been written down in ref. [17]. For us the 
following interactions leading to flavour- violating reaction v^N — ► r + ... are 
interesting: 



L' = (giLq c L iT 2 lL + gmu c R e R )S 1 + {h 2L u R l L + h 2R q L ir 2 e R )R 2 + h.c. (3) 
-m\StS 1 - m 2 2 R+R 2 

Here qi, II are the left-handed quark and lepton doublets, and e R , d R , u R are 
the right-handed charged leptons, down- and up-quarks, respectively; ip c = 
Cip T is a charge-conjugated fermion field. The subscripts L, R of the cou- 
pling constants denote the lepton chirality. The indices of the LQ's give the 
dimension of their SU(2) representation. Colour, weak isospin and genera- 
tion(flavour) indices have been suppressed. After integration over heavy lepto- 
quark fields one can obtain effective flavour-violating four-fermion interactions 
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W = ^(^r 2 / L )(g^r 2 / L ) + + 9 -l^{f L ir 2 l L e R u R ) (4) 

771 ^ 771 ^ 

. h 2L h 2R i - ■ \ , 
+... H 2 —(u R l L e R iT 2 qL) + ••• 

Especially interesting are couplings with the first and second quark genera- 
tions: 



T i 9lL,igiL, 2 /_ c _ c w-jc /r\ 

L 4AF = 2 ( S L^L - C L [i L ){v T d L - (5) 

7YI y 



- C \ , 9lL,l9lR,2 ,-c - c -c - c\ , 

T LU L ) + ~ 2 [s L V flL T R C R - C L fi L T R C R ) + 

771 ^ 

9lL,l9lR, 2 i T C - c -c - c\ , h 2L 2 h 2R i . 
2 [d L v ilL T R c R - u L n L T R c R ) H - 2 [T R d L c R v^ L - C R H L T R U L ) 

771 ^ /772 

^2L 1 ^2R 2 / — — — — \ 

H ' 2 ' [UrV^lTrSl - u r ii l t r c l ) + h.c. 

m 2 



For the couplings of Eq.(5) bounds resulting from flavour- violating r-lepton 
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decays bounds are absent. So, the CHORUS and NOMAD data are in fact the 
single ones that allow to extract the most stringent bounds on couplings of 
four-fermion interactions of Eq.(5). In quark-parton model the cross section of 
the inclusive u^N — ► /i + ... scattering on the isoscalar nucleon is determined 
by the formula Q 



a(z/ M iV -»• fj, + ..) = ~(<t{u,jP ->(! + ..) + oiv^n -»• /i + ...)) = (6) 

= ^((Q + £) + -(<? -£))- — -0-22 

Z7T O 7T 

where [16] 

Q = U + D + S= [ [u(x) + d(x) + s(x)]dx = f q(x)dx « 0.44 (7) 

JO JO 

and G is the Fermi constant. 



Q=U+D+S= f 1 \u(x) + d(x) + s(x)]dx = (8) 
Jo 

f q(x)dx « 0.07 

JO 

[/ « 0.03 (9) 

0.01 (10) 

Ufa 0.28 (11) 

73^0.15 (12) 

Using the four-fermion flavour- violating interaction of Eq.(5) and quark-parton 
model one can find that 



a(^N - r + ...) = ^[^^(35 + i(5 + ff)) (13) 

+ ^2|^ + 9 J L 2 mRa { U + D ) + 

mf m{ 2 

^^&U + -D)) + ^4^(5 + + 73))] 



Here we use notations of ref. [16] 
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From formulae of Eqs.(6-13) we find bounds 



9lL,29lL,2 



2 

mf 



<0A-G 



(14) 



9lL,29lR,2 



nil 



< 0.16 • G 



(15) 



9lL,29lR,l 



< 0.08 • G 



(16) 



hl^L < o.03 • G 



ml 



(17) 



h 2 L,ih 



2R,2 



ml 



< 0.16 ■ G 



Note that in supersymmetric models with R-parity violation [16] one of the 
possible R-violating terms in the superpotential has the form 

W = g 1L l^T2qLd c R (19) 

After integration over right-handed down squarks we find in particular the 
effective interaction 

L AF = 9lL ' 129 2 1L > 23 (-u c L VL + JL"hl) ■ {-clr L + r L v T y + ... (20) 

which in fact coincides with the first term of the interaction (5). Here m sq is 
the mass of right-handed down-squark. So using the previous results we find 
that 

glL.23glL.12 < Q Q8 . q (21) 



m 



l 



For 5 , iL,23fl l iL,i2 = 1 we find that m sq > 1000 GeV. 

One can extract bounds on effective coupling constants of the general four- 
fermion interaction describing the transition u^N — ► r + ... Q 

3 Here we consider only interactions of the first and second quark generations with 
and r 
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La = -^[G ljU u R iy^T R d L + G hl2 u R u^r R s L + G 1:21 c R u^T R d L + G 1:22 c r u i ,t r s^) 

G 2 , 22 s c L u^f L c c L + +G 2t21 s c L u^f c L u c L + G 2>12 d c L ^f L c c L + G 2M d c L u^f L u c L + 
G 3 ,ud c L v»f R u R + G^ 12 d c L v^f R c c R + G^ 21 s c L u^f R u c R + G^ 22 s c L u^f R c c R + /i.e.] 

In quark-parton model for the interaction of Eq.(22) the z/ M — > r inclusive 
cross secrion is determined by the formula 

- r + ...) = yf^lnt^ + + ^ + !>)) + Gj^U + D)) + (23) 
Gl 12 {\{U + D) + S) + Gl 22 (S) + (G| n + Gl n )^(U + D) + \{U + D)) + 
(C| 21 + G 2 21 )(3S + \{U + £>)) + {G\ 12 + G 2 12 )(^(C/ + £>)) + (G| 22 + G 2 , 22 )(3S)] 

Using NOMAD bound of Eq. (3) and Eq. (6) for inclusive — > /i cross section 
we find bounds on parameters of flavour-violating four-fermion interaction of 
Eq.(24) 

(Gi,n, G 2tll , G 3) n) < 3.3 • 0.8 • 10- 2 G (24) 

Gi,2i < 2.8 • 10~ 2 G (25) 

Gi,i2 < 5.5 • 10 _2 G (26) 

Gi i22 < 13.0 • 1(T 2 G (27) 

(^2,12,^3,12) <2.6-10- 2 G (28) 

(^2,12,^3,12) < 7.5- 10- 2 G (29) 

(^2,22,^3,22) < 7.5- 10- 2 G (30) 

It should be noted that in the search for /i-t conversion we probe in particular 
four-fermion flavour changing interactions. As it has been mentioned in the 
introduction bounds on the flavour changing interactions 

AL = jirlduc + G 2 bd + G 3 sb] (31) 

are not very stringent due to trivial kinematical reasons. For, instance r lepton 
cannot decay to D-meson and muon, since m T < m M + mc). So, bounds from 
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r lepton decays do not lead to bounds on Gj. Due to the same reason bounds 
from rare D- and B- mesons are also not very stringent. The study of ji-T 
conversion will allow to obtain stringent bounds on coupling constants G{. 



3 Experimental search for the jj, — > r conversion. 

Here we describe the search for flavour- violating interactions of the \ifqq type 
by using the muon to tau-lepton conversion on nucleons at high energies. Note 
that the four-fermion interaction \ifcu does not contribute to r-decay, since 
.D-meson mass is bigger than r-lepton mass, so for such interaction the flavour- 
violating bounds resulting from r-decay modes are not applicable, hence the 
search for jj, — r conversion is in fact the single way to probe such interactions. 

As an example, we simulated the experiment on searching for jj, — > r conver- 
sion performed with the NOMAD detector [20] at a high purity sign-selected 
20 GeV muon beam. The experiment is based on searching for single muon 
events from the primary muon interactions with a target (see Figure 1) with 
considerable energy losses which are not detected. The simulation were partly 
based on the Monte Carlo program used at NOMAD for the standard neutrino 
interactions [21]. 

The NOMAD detector, see Figure 1, consists of a number of sub-detectors 
most of which are located inside a 0.4 T dipole magnet: drift chambers (DC) 
with an average density of 0.1 g/cm 3 and a total thickness of about one ra- 
diation length (~ 1.0X ) followed by a transition radiation detector and a 
preshower detector not shown in Figure 1, and a lead-glass electromagnetic 
calorimeter (ECAL). A hadron calorimeter (HCAL) and muon stations are 
located just after the magnet coils. For simplicity, we have assumed that the 
primary muon interactions occur in a fully active dense target which is a block 
of lead glass, so that the energy losses in the target are measured. 

The sensitivity to fi — > r conversion and background level were studied for 
events with a simple topology of the final state. Namely, we consider quasi- 
elastic (QE) fi — > r conversion fi + N ^ t + N at a single nucleon N with 
the subsequent r decay in the target. For deep inelastic events the level of 
background was found to be considerably higher. 

Consider the r — > fiuu decay of r's produced in the /i — > r conversion. Since r- 
leptons are very short lived particles - even at energy of ~ 100 GeV the average 
decay length of r lepton is of the order of a millimetre - practically they are not 
detectable. The experimental signature of their presence in the beam would 
be a signal of fractional "disappearance" of primary muon beam energy in the 
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Magnet 




Fig. 1. Schematic illustration of the "missing energy" type experiment on production 
of t via n — > r conversion in the active target and its detection via r — > [ivv 
decay. The experimental signature of the fi — > r conversion is a single muon in 
the final state with a catastrophic energy loss in the target. The muon momentum 
is measured by the drift chamber (DC) spectrometer. The muon is accompanied 
by no significant activity in the electromagnetic calorimeter (ECAL) and hadronic 
calorimeter (HCAL). 

detector.^ Indeed, two neutrinos from the r — > \wv decay would penetrate 
any type of calorimeter without significant attenuation and cannot be observed 
effectively in the detector via its interactions. Hence, the only visible energy is 
the energy associated with the decay muons. The experimental signature for 
the fi — ► t conversion in the detector is a single muon of the same sign and 
with energy less than that from the primary muon beam, Figure 1. 

The energy distribution of muons from r —> \wv decays of r's produced in 
the fi — > r conversion is shown in Figure 2. The total efficiency of single muon 
detection is about 11% and dominated by the branching ration of tau to muon 
decay mode of 17 % and by the muon momentum cuts. 

There are several potential sources of the background for the experiment con- 
sidered: 

• low energy tail in the muon beam energy distribution, e.g. due presence of 
pions in the beam and their decays in flight; 

• catastrophic muon energy loss in the target with poorly detected products 
of the reaction, e.g. deep inelastic muon scattering accompanied by neutrons 

4 Another possible way of r identification is based on the algorithms developed by 
NOMAD for searches for r's produced through Va — > v T oscillations [21]. 
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Fig. 2. Simulated muon energy distribution from r — > /iw decay 0/ r 's produced in 
fi — > t conversion at energy = 20 Ge V. 

or nuclei break-up; 

• the detector is not hermetic in the forward direction; 

• muon missidentification. 

The largest contribution to the background is expected from muon inelastic 
photonuclear interactions in the target yielding a muon and neutral penetrat- 
ing particles in the final state (e.g. neutrons, K^, ■■) which are not detected, e.g 
due to the fact that the detector is not fully hermetic, thus making hermiticity 
to be the crucial parameter. The effect of leakages depends on beam energy. 
According to the simulation, the probability for secondaries to be mismea- 
sured decreases with the muon energy increasing. Additional way to improve 
the sensitivity is to use active target with a low threshold and detectors with 
rather good energy resolution capability. At higher energies when the r-decay 
length L T >1 mm, the active target could be an instrumented target, e.g. 
silicon micro-strip detector, to measure a non zero impact parameter of the 
secondary tracks from r decays with respect to the main vertex [23]. This 
technique being combined with the "missing energy" measurements will dra- 
matically improve the sensitivity of the /i — > r conversion search. Note that 
the target could be the existing STAR detector (instrumented Silicon Target) 
which was built and installed in the NOMAD detector [24]. 

The distribution of the reconstructed energy for muon-nucleon inelastic reac- 
tions is shown in Figure 3. The muon produced in the \x — > r — > \wv chain 
is defined as a single track accompanied by no activity in others subdetectors 
by using the following selection criteria: 
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Fig. 3. Distribution of total reconstructed energy in the detector muon-nucleon in- 
elastic interactions. 

• energy deposition in the target is consistent with the one deposited by a 
going through muon; 

• single DC track with momentum 3 < P < 10 GeV/c matched to the single 
track in the MUON stations; 

• no 7's in the ECAL with energy E 7 > 0.4 GeV; 

• no HCAL activity: total HCAL energy < 1.4 GeV. This cut serves as an 
HCAL veto; 

The probability P(/i — > r) for \i — > r conversion and the number N(fi — > r — > 
\wv) of signal events in the detector are related by 

a(uN -> rA) , , 

a in {fiN -> /i...) 

and 

A(/i -> r -»• ^i/i/) = A^ ■ cr(//A -> rA) • BR(r -> /w) ■ e ■ p ■ L ■ N A (33) 

where A M is the total number of muons, cx iri (/iA — > //...) ~ 3/i6 for the energy 
range discussed [22] and a(/j,N — > rA) are the muon-nucleon inelastic cross 
section and cross section of /i — > r conversion per nucleon, respectively^], 
BR(t — > yuz/z/) is the branching ratio for r — > /iz/z/ decay mode, e is overall 



5 For simplicity we neglect the difference between cross sections on proton and 
neutron, shadow effect, etc.. 
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Total number of muons 



Fig. 4. Simulation: 90 % CL limit on probability of the fi — > r conversion in the 
target calculated as a function of total number of muons passing through the target. 
Solid line is a polynomial fit to the Monte Carlo points, dashed line is its extrapo- 
lation to the region of higher number of muons. 

detection efficiency, p, L are target related factors and Na is the Avogadro 
number. 

Figure 4 shows estimated the 90% CL upper limit for the probability of fi — ► r 
conversion calculated based on the expected number of background events as a 
function of the total number of muons passing through the target. The study 
of the background in the region of high sensitivity — > r) < 10 -6 ) re- 
quires simulation of a very large sample of events resulting in a prohibitively 
large amount of computer time. The maximum number of simulated events in 
our sample was about 10 6 events. Assumption that the extrapolation shown 
in Figure 4 is valid down to P(fi — > r) ~ 10~ 9 is probably too optimistic, 
however seems has some promise. More detail study of the simulated back- 
ground shape is needed with significantly increased Monte Carlo statistics. In 
addition, the question how reliable are Monte Carlo predictions in the region 
of P(fx — > r) < 10~ 9 has also to be studied. 



We note that direct experimental searches for a signal from the [i — ► r con- 
version can be performed, for example, in the framework of the COMPASS 
experiment at CERN [25]. For muon beam energy region >30 GeV used 
in COMPASS the advantage is that the muon-nucleon inelastic cross section 
increases and the sensitivity to the fj, — ► r conversion rate might be expected 
to be of the order of magnitude high than shown in Figure 4. The search for 
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flavour changing reaction fj,N — > r + ... will allow to test flavour changing 
interactions at the level 10~ 2 G in terms of the corresponding Fermi coupling. 



4 Conclusion 

In this note, we derive phenomenological bounds on the strength of fv^qq four 
fermion flavour changing interactions. Such interactions can arise in models 
with leptoquarks, R-parity violating supersymmetry and additional flavour 
changing Higgs boson interactions. 

The bounds are extracted from the combined result of NOMAD and CHORUS 
experiments on search for z/ M — > v T oscillations, which constraints the rate of 
the r-lepton production in u^N scattering. Some of these bounds supersede 
limits from rare decays. 

We also propose to probe flavour changing interactions by searching /j, — ► r 
conversion on nucleons in a "missing energy" type experiment. This experi- 
ment could be performed at present accelerator, e.g. at high energy at CERN 
in the framework of the COMPASS experiment [25] or at the future neutrino 
factory. The instrumented target could be the STAR detector which has been 
installed and successfully tested in the NOMAD detector [24]. The estimate 
shows that one might expect a sensitivity to the probability of fi — > r conver- 
sion on nucleons to be of the order of 10~ 9 or even better. 
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